arXiv:nucl-th/0503072v2 29 Jun 2005 


The generalised relativistic Lindhard functions 


M. B. Barbaro^, R. CennB and M. R. Quaglia^ 

^ Dipartimento di Fisica Teorica — Universita di Torino 
Istituto Nazionale di Fisica Nucleare — Sez. di Torino 

Torino —Italy 

^ Dipartimento di Fisica — Universita di Genova 


Istituto Nazionale di Fisica Nucleare — Sez. di Genova 

Genova — Italy 


Abstract 


We present here analytic expressions for the generalised Lindhard 
function, also referred to as Fermi Gas polarisation propagator, in a 
relativistic kinematic framework and in the presence of various reso¬ 
nances and vertices. Particular attention is payed to its real part, since 
it gives rise to substantial difficulties in the definition of the currents 
entering the dynamics. 

1 Introduction 

The linear response theory, whose key ingredient is the Lindhard function 
(LF), is usually formulated in many-body frameworks like RPA, Landau 
quasi-particle theory and the like (see, e.g., Ref. [T]). Noteworthy, these 
simple approaches are often able to describe successfully a quite involved 
physics. 

The LF j2I, as originally dehned, is just the particle-hole polarisation 
propagator for a non-relativistic free Fermi gas (FFG) of electrons, and reads 
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is the Green’s function of the free electron. The factor 2 in front of (CD comes 
from the spin traces and is replaced by a factor 4 in nuclear physics (spin 
plus isospin). 

The explicit form of CD is known since 1954 j2] and, more recently w, 
has been expressed according to 


n°(go, |q|) = 2m [I[yi\q\, go + irj sgngo), |q|] + I[y{\q\, -Qo - ip sgngo), |q|]] , 

( 3 ) 

where 

i[yA^\]= L i {y^F + l{y^-kl)\og ^ , (4) 

(27r)^|q| ( 2 y + kp} 

n I X mqo |q| 

2/(|qUo) = -^ (5) 

| 4 | ^ 

being the West’s scaling variable [S]. The analytic extension when the argu¬ 
ments of the logarithms become negative is prescribed to be 


log(x ± i? 7 ) = log bl ± Z7r6'(—x) (6) 

and the imaginary part of the LF, namely the response to a scalar (-isoscalar) 
probe, is thus obtained. 

Actually many calculations of a fermionic relativistic response function 
are available — and hence a number of relativistic many-body computations 
have been performed. Still the need of analytic expressions for the relativistic 
generalisation of the LF, as an useful input for a variety of calculations, is felt: 
indeed the results presently available (see, e.g., Ref. [Hj, which completes a 
previous result [Zj, and Ref. jH]) mostly refer to the electroweak case. However 
we are also interested to the relativistic response to pions and p-mesons, 
which lead to different generalisations of the LF. Further, the excitation of 
nucleonic resonances need to be accounted for and, last but not least, in the 
quark-gluon plasma the case of massless particles (or one massive — an s or 
a c quark — and one massless) deserves some attention. 

In this paper we address a number of the above cases, showing that they 
can be handled algebraically in terms of only a few explicit functions. 
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The scope we pursue is to provide, as a hopefully useful tool for people 
involved in the field, a rather comprehensive description of the generalised 
LF in the relativistic case for nucleons and for 1/2 and 3/2 spin resonances. 
The limiting case quoted above will also be examined in detail. 

2 Setting the problem 

The channel dependence of the FFG polarisation propagator in a relativistic 
framework is more pronounced than in the non-relativistic case. We start 
with the generalised LF for a non interacting nucleonic system 

n“,(|q|, *) = -»Tr I -0j;S„{p + , (7) 

where 

s (p) = / ^(IpI ~ ^ - IpI) _ 1 

2Ep \po - Ep + ir] po - Ep - ip po + Ep - ip 
= (]/+m)6^(p) (8) 

= So{p) + ^^J^ 2i7r6{po - Ep)e{kF - |p|) 

is the nucleon [electron, quark...] propagator in the medium, with 

Ep = \/p'^ + 171^ , ( 9 ) 

the indices x and y label the incoming and outgoing channel (not necessarily 
coincident) and (9*, are some combinations of 7 matrices and momenta 
embodying the vertices characterising the channels. Moreover, Sq is the 
analogous of Sm in the vacuum. 

We have also introduced a “reduced” fermion propagator &m{p) spoiled 
of the Dirac matrix structure, its inverse reading 

T){p) = p^ — . (10) 

Although as given in Eq. (ED is ill-dehned since it is expressed by a 
divergent integral, we will not require renormalisability, since one may also 
be interested in effective theories. However, a regularisation procedure is 
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needed in order to cancel the divergences. In the case of o (one-loop level) 
the vacuum subtraction is sufficient p. 

Thus, dehning the following polynomial in the relativistic invariants 
p ■ q and 

fxyip, q) = Tr(/ (j[ + m)Oy , (11) 

n°, will read 

/ (J 4 l 

[Sm(p + q)&m(p) - S(,(p + «)So(p)| q) , 

( 12 ) 

where &o = {j/ + m) ^Sq. 

The frequency integration in m then reduces to the evaluation of the 
residua in the (say) lower half-plane, because along the half-circle at inhnity 
the in the denominators become irrelevant and the integrand vanishes. 
Thus the regularised 11°^ is given by 


n°j/(|qUo) = 


d^p 


X 


(27r)^ AEpEpj^q 

p -h q| - kF)9{kF - IpI) 


do Ep+q T Ep -|- ip 

6{kF - Ip + q|)6'(|p| - kF) 




do ^p+q T P^p 

e{kF - Ip + q|) 


ip 

-JxyiP^d) L„ = 


PO=Ep+q-qo 

fxy{p,d) Ip^^Ep 


go - Ep+q - Ep +ip 

_ - IpI) _ f („„) 

QO + -^p-\-q + -^p ~ 

+ 6'(|p + q| - kF)6{\p\ - kF) 


P0=Ep+q-qQ 


'po=Ep 

fxy{p,d) 


'po=Ep+q-qo 


-fxyiP^d) 


Po=Ep 


do ~ Ep+q + E. 


(13) 


(note that the term in the last line is never singular). 

Now in each denominator the factor ±ip can be replaced by +ip sgn(go). 
In fact the denominators in the hrst and third term can only vanish when 
go > 0 so here the replacement ip —>■ ip sgn(go) is immaterial, while the 
second and fourth ones can vanish for go < 0 so that the term —ip plays the 
same role of ip sgn(go). 
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Since, as we shall see, the explicit form of f^y is inessential for the fol¬ 
lowing discussion, we consider the case fxy{p,<l) = 1- Then, with some 
manipulations of the 6 functions and a change of variable, Eq. m simplihes 
to 



1 1 



(14) 


Eq. (fTT|l displays a great advantage from a practical point of view since 

1 . each term contains only one 6 function, hence the analytic calculation 
of the integrals is simplihed; 

2 . n° can be evaluated in a region where it is real and then its imaginary 
part follows by analytic extension by suitably approaching the real axis 
in the complex plane of go; 

3. it is manifestly even in go. 

The same procedure led to the form © for the non-relativistic LF. 

Now consider a N* excitation of mass M. In this case the polarisation 
propagator at the lowest order is built up by two non-coincident Feynman di¬ 
agrams, namely (the star will always denote quantities involving a resonance 
or a resonance-hole pair) 



(15) 


where 


6 *(p) = [p - + it]] ‘ 


(16) 


For later purposes we also introduce the inverse of ©*, namely 

S}*(p) = . 


(17) 
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The explicit form of will be specified later. 

Clearly, each term in m contains only one 6*-function because there is 
no Pauli blocking. Convergence is ensured by vacuum subtraction and, at 
variance of the nucleon-hole case, Eq. (jHl) (third line) tells us that po can be 
replaced everywhere by Ep, since So never contributes. 


3 Structure of the Lindhard functions 

In this Section we set up the general structure of the relativistic LFs, which 
will be later evaluated in some specihc cases. 

Before presenting the detailed calculation, we observe that the Lorentz 
covariance is broken by the presence of an inhnite medium like the FFG, 
since this naturally selects a privileged frame of reference, namely the one in 
which the FFG is at rest. Indeed in this system the nuclear matter has zero 
momentum, while any boost, no matter how small the velocity is, generates 
a state with inhnite momentum. 

If we instead consider a system with mass M and hnite momentum p, then 
any response function / to a probe carrying a four-momentum = {go, q) 
can only depend upon Lorentz scalars, namely / = f{p^ = M‘^,p ■ q,q‘^)] 
however in the M —>■ cx) limit p - q = qo \/p‘^ + ~ P ' q so that the 

Lorentz covariance is broken since / will depend upon qo and q separately. 


3.1 The ingredients 

Having clarihed the functional dependence of the LFs, we now introduce the 
ingredients needed to their evaluation. We dehne the functions 




d^p ejkp - IpI) 

2Ep {p + qY — -|- ip 


to be computed in the next Section. The quantities of direct physical interest 
are the even and odd parts (in qo) of , namely 


|q|,g„) = U, 


_ jj*{n] 


rel 


(M; |q|,go) ± |q|. - 1 o) ■ 


(19) 


We shall consider in the following a variety of functions f^y, each one giving 
rise to its own LF: remarkably these all are expressed in terms of few basic 
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cases. Note that when f^y = 1 then Eq. m becomes 


n°(|qUo) = |q|,go + * sgngo) • (20) 


The above functions display a complex analytic structure (logarithmic cuts) 
and have a well dehned imaginary part and a well dehned asymptotic be¬ 
haviour (in qo), namely ~ this means that the real part of 
hence of the LF, can be univoquely recovered from its imaginary part via 
dispersion relations. 

However in general f^y is a polynomial in gg. Each term of this polynomial 
generates contributions to n° with the same imaginary part (up to trivial 
coefficients), but with different asymptotic behaviour, so that the evaluation 
of the real part will require subtracted dispersion relations. The subtracted 
parts will be called contact terms. These can be expressed in terms of the 
fc^T’-dependent function 


= -i 


(^ 


<5m(p) 


renormalisation 


d^p 9{kF — IpI) 


Stt^ 


-' J (27r)3 

kpEp — rn? log 


2Ep 

kp -|- Ep 


m 


( 21 ) 


where Ep = \/kp + m? and, more generally, 

^\n] [ d?P ^(^f-IpI)™-! p /3 1 n_5_ kl\ 

( 22 ) 

where 2 E 1 is an hypergeometric function. 

Note that contact terms are also indirectly related to the renormalizability 
of a theory, because their presence alter the power counting in a bosonic 
loop entering in an RPA-dressed bosonic propagator closed on itself or on a 
fermionic line. 

Now we are in a position to deal with the general structure of the LFs. 
Since we shall consider (pseudo-)scalar and (pseudo-)vector couplings, our 
LFs will carry 0, 1 or 2 vector indices only. Tensor couplings could bring into 
play other functions and but they seem not to be, at present, of 

physical interest. 
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3.2 Scalar case: 0-index functions 


Here and O'^ have no vector strnctnre, hence f*y = fxy{p^,P ■ in 

general a polynomial, is a Lorentz scalar. Replacing then with '£){p) + m? 
and nsing the identity 

p-q = \[{p + qf - - \p(f 

1 ^ (23) 

= 2®*(P + ^) “ - 2^^^ ’ 

we can rewrite f*y in the form 

f:, = {®’(P + «)}"■ {®(P)}” • (24) 


being A^{q^) and A^„(g^)’s Lorentz scalars (hence the snperscript S). In 
(OH|) we have introduced, as in HU, the dimensionless qnantity 


P 



(25) 


In the case M = m we have f*y —>• f^y and p = 1. Fnrther, in Eq. the 
snmmed indices m, n mnst satisfy m + n > 1 . 

If we consider the resonance-hole case, we have to insert (121 into dn. 
Then the first term on the r.h.s. of Eq. (j21 yields Aly{q^)T*^\M] |q|,go) 
and the second generates the contact terms, which in the present case can 
be expressed in terms of the following fnnctions 


rA”’”'(|q|.9o) = -i j (04 {S-(p + «)}’"-' {B(p)r‘ 

=F*/( 0 i{®*(p- 9 )r-‘{B(p)r‘ . (26) 

Then, becanse of the vacnnm snbtraction, only the n = 0 term snrvives in 
(121 (otherwise any dependence npon kp is lost) and of conrse it mnst be 
m > 1. 

The 3/2-spin resonance generates an additional complication dne to the 
possible presence of projection operators, as discnssed in Sec. |H] (see Eq. (IM|l 
for details): these will reqnire the addition of a fnnction B^{q‘^), whose role 





will be clarified later. In conclusion, the most general 0 -index LF has the 
structure 

|q|,go) = |qUo) = |qUo) 

+ = 0; |qUo) + A^o(9^)^s!r’°’(|qUo) 

m>l 

(27) 

with A^y{q‘^), B^y{q^) and A^Q(g^) to be specified according to the problem 
one deals with (however B^y{q^) = 0 for spin-1/2 particles). 

The functions Tg|("’°^(|q|, go), linked to the functions T of Eq. ( 1 ^ . are 
not Lorentz invariant. Those entering our calculations are explicitly given in 
Appendix 

The nucleon-hole case is more involved. Eq. (El still holds valid, provided 
T)* = D, but the subtraction scheme will be different, since D{p + q) also 
depends upon kp- Thus the contact terms will also be different because both 
the cases m = 0, n 7^ 0 and m 7^ 0, n = 0 contribute in this instance after 
the vacuum subtraction. 

3.3 Vector case: 1-index functions 

Vector-like LFs can only arise through the combination of a scalar and a 
vector vertex. Lorentz invariance would forbid such transitions, because 
scalar and vectors belong to different representations of the Lorentz group 
but, since the infinite nuclear medium violates covariance due to the presence 
of the 6*-functions, these terms may occur. Hence in the nuclear medium a 
vector meson (the u, for instance) can be converted into a a. 

By covariance the functions f^y must have the structure 

/4 = + ( 3 q^ , ( 28 ) 

where the transverse momentum 

( 29 ) 

qz 

has been introduced {t ■ q = 0 ). The second term on the r.h.s. of Eq. (OH|) is 
immediately handled, since the vector q^ factors out of the integral and the 
scheme of the previous subsection applies, but it gives rise to a H^y which 
is not gauge invariant. 
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Instead, the vector-like LF generated by the hrst term of f^y (to be called 
cons^ obeys the conservation law = 0 and can be cast into 

the form 


TT*(l)/i cons _ TT*(1)0 cons^/i. 


(30) 


where we have introduced the four-component object (not a vector) 


= 



(31) 


Hence it is sufficient to compute the 0 component of H 


xy 


only. Dehning 


(M; |q|,go) 




(27r 


{e*{p + q)em{p)±&*{p-q)&Up)} (32) 


and 


T'vr’^klqUo) =-* J ^±{qo^-qo) (33) 

the 0 component of the vector-like LF, using again (OH|l . takes the form 

= H^^(g2)0f (M; |q|, go) + Bl{q^)£lf {M = 0; |q|, go) 

+go) 1 

m 

(34) 


where we have accounted for the projection operators (inn|) and the parity is 
not specihed. Again the functions A, B and A have to be specihed according 
to the problem. 

The new function £l^ is expressible in terms of the T*!”! as follows 


nf(M;|q|,go) 


T±W(M; |q|,go) + ^pgoT^^°^(M; |q|,go) 


0 

^tj;[0] ) (35) 


while the relevant to us are listed in Appendix^ In conclusion the 

general structure of the 1-index LF reads 

n:«'^(M;|q|,go) =n:J'+(M;|q|,go)^^ + g'^n:J_(M;|q|,go) . (36) 
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Again the nucleonic case is more tricky, because we must first replace 
with 



( 37 ) 


next express p ■ q according to (I2d|l and hnally use directly the expression 
ca. Following exactly the same path leading to Eq. (HI we obtain the 
limiting case (M = m) of Eq. (jHSl)- 

3.4 Tensor case: 2-indices functions 

The 2-indices functions, which require two vector-type vertices, can be split 
into a symmetric and antisymmetric part that need to be studied separately. 

3.4.1 Symmetric case 

In the symmetric case, since f^y is a true tensor, Lorentz covariance imposes 
the structure 


= Oi ilV' - rit) + + as(tv + qV) + (38) 


the Oj being Lorentz invariants. 

The hrst two terms of correspond to a conserved current. Let us 
denote the associated LF by Llxy^lymm" '^(M; |q|, go), which can be split into a 
longitudinal (LI*^) and a transverse (n*^ polarisation propagators, dehned 
according to 


n*^(M;|q|,go) 

n*^(M;|q|,go) 


tt*(2)00 cons 
symm 


[M] |q|,go) , 


(39a) 

(39b) 



and, in a compact notation. 



(M; |q|,go) 



( 40 ) 
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The first term in the r.h.s. of (EHI) is easily handled, since it reduces to 
the 0-index case. The second instead requires the introduction of two new 
quantities with the associated contact terms. We thus define 


nf(M;|qUo) 


d^p {ty e{kF-\p\) 

{2'kY 2Ep {p + qY — + ip 


= \q\,qo) + qopT*^^\M] |q|,go) 

1 f 

+ |qUo) - < 



± (g <-^ -q) 

PQ = Ep 


(41) 


and 

nf(M;|qUo) 


1 

r d^p 

IpI 

2 

|q| 

1^ - (p ■ q)^ - 

IpI) 

|q| 

N 

(27r)3 




2Ep {p + g)2 — -I- ip 


Po—Ep 


± (g ^^ -g) 

_ Q \ 'y>*[ 2 ] 


|q| 


|q|,go) + gopT*'^^(M; |q|,go) 


^ q^P 


+ I '^q^P'^ - 1 |q|,go) + < 


2 |qr 


X[o] 


together with the contact terms 

rh”’"'(|qUo) =-« j {s-(p + 9)r-‘{®(p)r-'±(«o 


(42) 


■9o) 

(43) 


and 


I?L"'"'(lqU„) 


|q| 


<iV 


Op 


2 | |2 

|q| 


IP 


q)1 


X 


(44) 
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{^*{p + q)r-^{^{p)r-^ ± (go - -go) • 


Thus, applying (OH|l and (EH), we obtain 

|q|,go) = -|q|^n;J(M; |q|, go) + A^J(g^)n;^(M; |q|,go) (45) 

+ B^y (g^)n^(M = 0; |q|, go) + ^ A^^(g2)T;|("’°^(|q|, go) 

m 

n*J(M; |q|,go) = -2q‘^n*J^{M; |q|, go) + v4^J(g^)l}(f (M; |q|,go) (46) 

+ B^'^{q^)Q*^{M = 0; |q|,go) + ^ A^'(;(g2)Tyf’°'(|q|, go) 

m 

(note that the same coefficients and A^^ enter in both 11*^ and 

n*^). The above relations give the structure of the longitudinal and trans¬ 
verse LFs and thus fully describe through (lin|) . Finally, the re¬ 

maining terms of dSHl) can be reduced to simpler cases and one gets the hnal 
result 

|q|,go) = |q|,go) 

+ (g^Tt^ + g^Tt^)n:^_(M; |q|, go) (47) 

+ gVriS+(M;|q|,go) . 

The last term on the r.h.s. of (03), namely corresponds to a LF with 

the same structure of but with different ingredients: these will be called 
and A;^o, respectively, to avoid confusion. 

Again the nucleon-hole case must be handled separately, by using the 
analogous of Eq. 03- Clearly now, at variance of the resonance case, only 
the terms exist, and a straightforward calculation shows that Eqs. 0T| 

and (1^3 still hold valid. 

3.4.2 Antisymmetric case 

An antisymmetric tensor should have the form 

= '>1 (*'‘9’' - 9'“*“) + (48) 

and the general structure of ^*}y'^a^tisyrara accordingly be 

n:™i.,„..,(V;|qUo) = {g“<n‘'-q‘'W)n:';jM-,\q\,g„) 
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+ • ( 49 ) 

Again the function n*y_ entering the above has the same structure given by 
Eq. (inH), but with the functions and A^q replaced by and 


4 Analytic evaluation of |q|,Q'o 


In this section we explicitly compute the function f/*]; dehned by Eq. (IT^ 
(the other two functions and Ulf} can be obtained along the same path 
and the results are reported in Appendix lB|l . 

Assuming here 9= go 7 ^ 0 we get 


f d^p 1 6{kF — IpI) 

J (27r)3 2Ep 2Epqo - 2p ■ q + g^p 

kp 

1 f gV+ 2p|q| + 2qoEp 

levr^lql y Ep g2p - 2p|q| + 2goEp ■ 
0 


(50) 


Note that the dependence upon M is fully embodied in the inelasticity pa¬ 
rameter p. Integration by parts yields 


qo] ^ 1 jp , gV + 2fcjr|q| + 2Ej.go 

Ibvr^lql ^ ^ g^p — 2/cir|q|-|-2Ep’go 

kp o o (51) 

_L ^ P^P ^ _ 

47r2 J ^ {q^p - 2p|q| 2Epgo) (gV 2p|q| -h 2Epgo) 

0 

The integrand in m displays four poles, located at p = ±p^. Dehning 

Q± = (M±m)^ —g^ (52) 


and 

A* 



VQ+Q- 


\/[{M — my — g2] [(M -|- m)^ — g^ 

g2 

( 53 ) 
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(with the chosen sign A* is positive in the space-like region) we obtain for 
the poles the expression 


A = Y'’ I 

The are the branch points dehning the region where, for real (ED 
develops an imaginary part. In particular the lowest positive branch point 
is just the lowest possible longitudinal momentum for the occurrence of a 
resonance-hole pair and consequently coincides (up to a sign) with the y 
scaling variable. Indeed —y*_, taken at M = m, is just the y scaling variable 
for the relativistic Fermi gas uni. 

The integrand of (ED can be split according to 


+ ^ _ ^ _ 

327 r 2 |q| (p + y* p + yX P - y*- P - yX} 

I { R*_ RX R*. RX\ 

^ 327 r 2 |q| \p +y* p + yX p-y*-^ p-yX) ' 


where 


RX = 


|q| 


A* ± 


Qo 


-P 


( 66 ) 


has the property 


{RXf =m^ + {yXf (57) 

and is trivially linked to the scaling variable ^|J* used in Refs. dniiaiisiiiii 
by 

R*_ = m + (Ep — m)-^*^ . (58) 

We can now easily compute explicitly ED, getting 


rr*[ 0 ]/ II 6 Pi kp + Ep 

Urel iM; |q|,go) = -y^^log 


q^p 


(47r)^ m 

■ 2 /c^|q| -|- 2,EpqQ 


Ibvr^lql ^ q'^p + 2kp\q\ + 2EpqQ 

qoP ^ {kp - y*_Y{m^ + kpyX + EpRX){m^ - kpyX - EpRX', 
647r2|q| ^ {kp + yXYirn^ + kpy*_ + EpR*_){m? — kpy*_ — EpR*_] 

m\kp - yXf{kp + yX? 

[m^ - {kpy*_ + EpR*_Y] [m^ - {kpyX + EpRXY] 


( 59 ) 


A 
647r2 
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or, with some algebra, 


TT*MfT\j II ^ y-+y+^ kp + Ep R*_^_ — R*_ + 2Ep 

Urel i^l |qUo) = - .^_ 9 I_T log - + -X 


IbTT^Iql 


m 


647r2|q| 


{kp - y*_fim^ + kpyl + EpRl){m^ - kpy^ - EpR^) 
^ {kp + yX)^{m^ + kpy*_ + EpR*_){m? - kpyX - EpRX) 


RX + R*_ 
647r2|q| 


log 


m^{kp - y*f{kp + ylf 


- {kpy*_ + EpR*_y] [m^ - {kpy^ + EpR^^] ’ 

-1 


which depends, but for the overall factor |q|~ , only upon the two scaling 
variables yX- 

Finally, with i?* = (hence E^^^^ = \/+ {kp ± qY) and 

bringing back go to the real axis, Eq. reduces to the compact form 


P^rel IqIj^o) — ~ H-or,_9i_.i ^2 “ 


Ibvr^ 327r2|q| 
having dehned the logarithmic functions 
kp -|- Ep 


647r2 


(61) 


h = log- 


h = log 


h = log 


+ i'Kk'i 

if A* is real and 


m 


% + ^kp-q + ^ 


go + + Ep 


+ log 


go ~ ^kp-q + E- 


+ mk2 


{2kp + goA*)^ — 

q 


{2kp - goA*)2 - 

q 



go - El^^^ + Ep 

^ {pkp - EpA*f - 

+ log 


(62a) 

(62b) 


4|q|^m^ 


{pkp + EpA*f 


4|q|^m^ 


k = 2*arctanh£^ _ 


di 


(^2 

,2 r,j2 


(62c) 


( 63 ) 


- 2i7r6'(cii)|6'(go) + 6{-qo)0{M^ + - 2kp + jq]^) 

X sgn(go + \J+ m? - 2kp + |q|^) 
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+ 2in9 (c/ 2 ) |6'(go) + sgn | go - 

cii=44-g2|A*|2-|q|V 

c/2 = + 


'M^+\qf+ "" 


+ 2k‘i 


4m^|q| 


(64) 

(65) 


if A* = /|A*| is purely imaginary (i.e., for < <il < i.Q’+Y^ see Eq. (|F)H|l i. 

In this case the function is manifestly real. The two 6 functions force 
/a to be continuous. In Eqs. (El k2 and /ca are integer that can be fixed 
by analytic extension or by checking the integral in some suitable points. 
They will be specified in the next section and found to depend only upon the 
analytic structure of the logarithms and thus remain the same for the whole 
set of functions 


5 The response region 

As previously mentioned, the singularities of El, independent of fxy, fully 
determine the response regions for particle(resonance)-hole(antiparticle) ex¬ 
citations. Further, the values of the ki in El are also independent of f^y 
and will be determined in this section, where we consider a real go up to a 
vanishingly small imaginary part ±/p. 

To set the response regions first consider the branch points associated 
with the excitation of a resonance in the free space. They are located at 

Q± = \J |q|^ + (M ± m)2 , (66) 

which dehne the boundaries of the regions where the production of an anti- 
resonance-particle plus the emission (go < —Q+) or the absorption (go > Q+) 
of a probe or the excitation of a particle to a resonance < go < Q^) 

are allowed. Accordingly and are usually referred to as the threshold 
and pseudo-threshold, respectively. These singularities clearly stem from the 
presence of the Dirac sea. 

Next we discuss the domain where the response of the system to a probe 
accounts for the effect of the medium, hence the effect of the Fermi sea. This 
is fixed by the logarithmic singularities of U*^i, which turn out to be located 
at 

Qijp) = - E, , (67) 
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that, for p = kp, fix the boundaries (hence the label B) of the resonance-hole 
and antiresonance-hole regions. An easy check yields the ordering 

> {Q^-? > {Qir > {Q^-? > ( 68 ) 

(uniformly in p). Also, it is found that 

Q+±{p) > -|q| , Q-±{p) < -|q| Vp (69) 

and 

iQir > |q|' • (70) 

Note that six critical values of |q| exist, corresponding to the various 
intersections of the boundaries. We have 

{Q--f = {-Qlf for q = ^ ^ ^ kp , (71a) 

m 

{Q++)^ = {Q-f for q = q^^ = ———kp . (71b) 

m 

Concerning the sign of the singularities one hnds that < 0 and Q++ > 0 
for any value of kp. Instead, one has Q^_ > 0 if kp < kp with 


k^^ = VM2 - m2 , 


(72) 


while, when kp > kf, is negative in the interval q^' < q < qcV, being 


(3) 


,( 4 ) 


r(3) = 

cr 

kp — A/m2 — M2 + kp = kp — y 

J^F 

j - 

{KY , 

(73a) 

rW = 

cr 

kp + y/TTl? — M2 + kp = kp + y 

J^F - 

to 

(73b) 


This occurrence deserves a comment; in fact, in the limiting case M = m, 
kf = 0 and we are always in the second case, with qif = 0 and qit'^ = 2kp. 
In this region the response function at hxed |q| and as a function of qo has a 
discontinuity in the derivative. 

Finally other critical points arise in connection with the light cone. It is 
obvious that (Q±)^ > |q|^, hence we must only inquire about We find 




for q < qir^ 
for q > gcr^ 


( 74 ) 
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where two new critical points appear, namely 


q 


( 5 ) 


q 


( 6 ) 


{M^ - m‘^){EF - kp) 
2m? 

- m^){EF + kp) 
2m? 


(76) 

(77) 


Now we discuss the response regions in the (g + 0, ||q||) plane and £x the 
corresponding values of k 2 and k^ appearing in (P|). 

Consider hrst the negative time-like region. Here for go < —Q'+ the 
antiresonance-particle production is allowed in the vacuum, but is ruled out 
by the renormalisation, which subtracts out the vacuum effects. It is, in 
any case, Pauli-blocked by the Fermi sea in the region spanned by go = 
~ Ep for IpI < kp, namely 

max{(5:^_(p)} > go > . 


The boundaries of the permitted response region are then found to be 


Q-pikp) < go < 


Ql 

Q^-ikp) 


for q < 
for q > qil^ 


(78) 


The resonance-hole region instead corresponds to go > — |q| and lives 
partly in the space-like region and partly in the time-like one. The allowed 
values of go span the interval 

min{Q+_(p)} < go < max{Q++(p)} , 


again with |p| < kp, the associated boundaries being 


Qp-ikp) <qo < 


Q^- 

Q^+ikp) 


for g < qir^ 
for g > qir^ . 


(79) 


In order to £x k 2 and k^, which are integer and constant inside all the 
response regions, it is then sufficient to evaluate the imaginary part of U*^i 
in some particularly simple point of the response regions. 





Consider first the regions (5++ < qo < and Q^_^_ < go < —Q+- Here a 
convenient point is |q| = 0, where 

^ / (27^3 2E ~ P)'^(2^p?o + pql) 

^ |q|=0 

which is non-vanishing only in the regions quoted above, its value being 


'^reZ-T 


|A* 


IGtt 


|q|=0 


This result deserves a few comments. First, from the dehnition it follows 
that the imaginary part of must have the sign =p according to whether 
Qo ^ Qo ^ iP) in accord with the above outcome. Then, since A* is negative 
in the time-like region, we obtain 

A ;2 = 0 (80a) 

/c3 = =f 4 (80b) 

for (5++ < go < Q- or < Qo < -Q+- 

Next we consider the region Q^_ < go < Q++- Here, choosing a very 

small value for kp, the point go = \J -1- |q|^ — m surely lies in the desired 
region, and one gets 



kp 1 

j - m){Ep - m) - 2pqx) . 

0 ^ -1 


Being kp (and hence p) small, the argument of the ^-function vanishes at 

P 


X = 


2|q| 


A _1 

m 


which is less than one. Thus the angular integration becomes trivial and we 
get 




U*[Q]rel = T 


Ep — m 


167r|q| 


Since this result should be a combination of the coefficients of I 2 and /s then 


k2 = Tl 


(81a) 
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h = . (81b) 

Finally we consider the region Q^_^_ < Qq < Q^_. Here we choose 
go = —E* — m and, following the same steps as before, nnder the same 
assumptions, we hnd that Eqs. m hold again. 

The singularities, and the response regions, in the plane |q|,go are dis¬ 
played in Fig. ^ (left and right panels) for two different values of kp, one 
below (left) and one well above (right) the critical value of kp- 


6 The limiting cases 

Having determined the response regions and the integers ^2 and k^ we are 
now able to derive the response functions in the most general case. In this 
section we consider some specihc limiting cases. 

6.1 The case M = m 

Here 

Q- = |q| , (82) 

Q+ = \/ |q|^ (83) 

and the critical points qcP and gc?^ occur at 2kp and 0 respectively. Further¬ 
more, kp = 0 and thus Q^_{kp) = Ekp-q — Ep is always negative between 
q% = 0 and g^j, = 2kp. Finally the Q^^{kp) always live in the space-like 
region. 

The expressions for U*^i given in Sec. 01 are still valid, provided we set 
M = m, which implies A = A* = and p = 1. Concerning the 

imaginary parts, nothing changes in the particle-antiparticle domain, while 
in the particle-hole region, since q^^ vanishes, the response is conhned to the 
range 

Qp-{kp) < go < Q^j^ikp) 

and k 2 and k^ are given by (EH). 
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6.2 The case m = 0 

This case may correspond to the excitation of a light qnark to an s or c qnark 
in a Qnark Glnon Plasma. 

Here 

Ql = e; , (84) 

while and qir^ tend to inhnity: accordingly and Q^_ never coincide. 
Fnrthermore kp = M. Finally it is immediately seen that gc? —> oo, implying 
Q++ > |q| ykp, iff, while = M‘^ /[Akp) . Thns the response region is 
represented by the intervals 

Qp-{kp) < go < Q+p{kp) 

Q^p{kp) < go < Q^_{kp) , 

where Eq. m holds, and by 

Q++{kp) < go < Q- 
Q--{kp) < go < -Q+ , 


where instead dHOI) is valid. 

Concerning the LF, since now p = A* = 1—^ and y^ = ^ j (|q|± 

go), we end up with the expression 


U, 


rel lm =0 


{2kp + go)g^ - M^gp 
327r2|q|g2 

m 2 f. (M 2 - 


log 


m2 - 2kp{qQ + |q|) - g2 


M 2 


327r2g2 


log 


2kpqo - q 


2\2 


2kp{qo — I 
4fc|,|q| 


q|) 

2 


+ ivrfcj 


(M2 - g2)2 


+ inks 


(85) 


6.3 The case M = m = 0. 

Finally we consider the extreme situation where both masses vanish. Here 
the values coincide with the light cone, while IhII ~ kp 

and inside the response region only the case of Eq. m occurs. 
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The expression of U*^i further simplihes to 


u: 


2kF + qo 


rel \m=0 



+ kpiqo + |q|) 


+ 2kF{qQ - |q|) 
— {2kF + qoY 


+ i'nk2 


( 86 ) 


+ iTiks 


7 The spin 1/2 resonances 

In this Section we explore the LFs associated to the excitations of the nucleon, 
addressing hrst the simpler case of the spin 1/2 resonances (e.g. the Roper 
(jV*1440) resonance). For sake of simplicity we disregard isospin, which 
simply yields a numerical factor. 


7.1 The 0-index functions 

Here the vertices we consider, beyond the identity, carry some 7 -matrix struc¬ 
ture of the type ^ or ^ times, eventually, a 75 . Owing to the mass shell 
condition for the nucleon, ^ — m) + m can be replaced by m, since 

/ — m cancels with the nucleon propagator, leaving a /ci?-independent term 
subtracted out by the renormalisation. The nucleon-hole case requires a 
separate discussion. 

Similarly, ^ = (^-|- ^ — M) — — m) — (M — m) leaves us with the 

identity times M — m plus, however, a contact term, an occurrence reflecting 
our ignorance about the off-shell reaction mechanisms. Actually the (][ vertex 
is redundant as far as the imaginary part of the LFs, and hence the response 
functions, are concerned. The real parts instead are altered by an extra 
contact term that matters in the response when higher orders (say, a RPA 
series) are accounted for. 

In attempting to account for the different off-shell behaviour one meets 
a proliferation of complicated and mostly irrelevant terms. Thus here we 
conhne ourselves to consider only the vertices 0 = 1 (a-meson absorption), 
75 and ^75 (pion absorption within the pseudoscalar and the pseudovector 
coupling). This last, at variance of the pseudoscalar coupling, correctly de¬ 
scribes the 7 r° suppression in the photo-production process and respects the 
chiral limit, owing to a further contact term added to the pseudoscalar vertex. 
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Notice that the vertices containing a 75 are derived from the corresponding 
parity-conserving ones (up to, eventually, a sign) by replacing m with —m. 
For the 0-index LFs Eq. (EH) applies and we only need to specify and 
in the various cases. 

Considering Erst a scalar probe, the function fss reads 

fss{p,q) = Tr(]/-f </+ M){^+m) = 4:{p'^ + p ■ q + Mm) 

= 2^*{p + q)-2^{p) + 2Q+ , ^ ^ 

where use has been made of the identity (El in the second line. Thus 
we get = 2Q+ and Afg = 2. The term 2S)(p) is cancelled by the 

renormalisation when studying the resonance-hole case, but it survives in 
the nucleon-hole one. In conclusion 

n:f(M;|qUo) = 2T5t’°'+2Q+Tf^(M;|qUo) • (88) 

The other vertices are decoupled from the identity by parity conservation. 
For the pseudoscalar coupling ( 75 ) we hnd 


n;g)(M;|q|,go) = -n:f(M;|q|,go) 


m—m 


For the pseudovector coupling, introducing the notations 


(89) 


Mt = M + m 5m = M — m , (90) 

we hnd A^ = 2M|.Q_ and Afg = 2{5mMT — q'^p), Afg = 2. As already out¬ 
lined, in going from the pseudoscalar to the pseudovector LF the coefficient 
A^ is multiplied by the expected factor M|., but the contact term does not. 

The pseudovector coupling conserves the axial current (or alternatively 
the existence of the Goldstone boson). Covariance would entail flpy ~ 
but, since it is broken, the Goldstone theorem only requires 


lim U*py{M-, |q|, go) = 0 . (91) 

qo^O 

IqHO 


Now we observe that 


(M; 0, 0) = 2Uf^ (M; 0, 0) = 2 


d^p 1 6{kF — IpI) 

(27r)3 2Ep m^-M^ 


2 q;[o] 

M‘2 — 

(92) 
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and it is easily verified that the contact terms are tailored in such a way to 
exactly cancel in the limit g ^ 0 . 

Finally the mixed pseudoscalar-pseudovector LF also exists: it has = 
2MtQ- and Afg = 26m while the pseudovector-pseudoscalar function has 
the opposite sign. 

The nucleon-hole case must be considered aside because of the different 
structure of the contact terms. For example, the scalar-scalar response, owing 
to (jHZl) leads to the manifestly vanishing contact term 

It is found that the only existing contact term pertains to the pseudovector- 
pseudovector case and is given by 

7.2 The scalar-vector interference 

In inhnite nuclear matter a scalar probe can be transformed, through a reso¬ 
nance (nucleon)-hole propagator, into a vector one. Thus we shall consider, as 
before, the scalar-type vertices I and 75 . Concerning the vector-like vertices, 
24 independent currents exist, 12 of them parity conserving and 12 parity- 
violating (see, e.g., Ref. US ). However, currents embodying a (that can 
be extracted out of the integral) times a (pseudo-)scalar structure reduce 
to a 0-index LF, already handled in Sec. 17.IL times q^. Further, neglecting 
contact terms, ]/ and ^ are redundant and the Gordon identity 

^ + ^ 7 '^(j 6 - m) 

+ i(M + m)Y - (93) 

allows us to express the current in terms of the usual currents 7 ^ and 
(j^^qy only on the mass shell, while the off-mass-shell extension of the currents 
remains unpredictable. 

Disregarding the huge variety of contact terms and q^, only four inde¬ 
pendent currents survive, and they may be forced to be conserved by adding 
some suitable terms proportional to q^. They read 
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jS = (94b) 

is ,4 = if,275 • (94c) 

Furthermore is expressible in terms of the other three currents by exploit¬ 
ing the charge conjugation symmetry. 

Since all these current are conserved, only the hrst term of Eq. dSHD is 
required. The functions are listed, with a self-explanatory notation, 

in Tabled 



I - Jf 

I - 

■ a 

■U 

75 ^ J4 


AMx 

2 2 
—r 

m 

46m 

2 2 
—r 

m 


Table 1: The function A^y{q‘^) for the scalar-vector Lindhard function 


7.3 The 2-indices response 



3i 

32 

fi 

00 

0 

32 

0 




3i 

32 

jr 

2Q. 


q2 

2 — 
mq^ 

32 



9 ^ ~ 

mq^ 

2m2g2^" 


Table 2: The functions A^{q‘^) for 
the vector-vector parity conserv¬ 
ing currents 


Table 3: The functions A^'^{q‘^) 
for the vector-vector parity con¬ 
serving currents 


We consider now the vector-vector response and distinguish between three 
different sets of LFs. 

1. First we examine the parity conserving-parity conserving LFs (set up 
by the conserved currents and j^), which are symmetric tensors. 
Thus only the functions 11*'^ and 11*^ are required, that in turn need 
the knowledge of A^, A^'^ and of the corresponding A’s, according to 
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(gSl) and (Hn|l . The are summarised in Table H and the in El 
The contact terms Afg pertaining to are displayed in Table 01 The 
other contributions, namely the A^^(g^) are all vanishing. Instead a 
coefficient Afg survives for the case and the relation \ 20 U 2 j 2 ) — 
l/(2m^g^) holds. 



3i 

h 


2 

q2 

5m 

mq^ 

42 

6m 

q^p — MT6m 

mq^ 

2m^g^ 


Table 4: The contact terms Afo(q'^) 
for the vector-vector parity con¬ 
serving currents 



2i 

J 2 

23 

Ai 

Mj' 

2i 

m 


5m 

97 

,MT5m 

^ t 

m 

^ 2 


Table 5: The functions A*^{q‘^) for 
the vector-vector parity conserving- 
parity violating currents 


2. Next we consider parity conserving-parity violating LFs (currents 
and j 2 at the incoming vertex, jg and at the outgoing one). Here 
the tensors are antisymmetric and Eq. (I49|l applies with the second 
term only, namely 


with dehned by Eq. (EH). The required functions A*^[q^) are 

displayed in Table El No contact term exists. The case would 
provide a non-vanishing function A]^, but actually = 0. Finally if 
we reverse the vertices the following relation holds: 




m+lTT(2)/ii^ 


^_l)m+in 


1 ^ 
JmJn 


m = 1, 2 , n = 3,4 . 


(95) 


3. The parity violating-parity violating LFs are derived from the parity 
conserving-parity conserving ones through the simple relations 


n* 

J3,J3,4 

n* 

h/4J3.4 


n 


■Jl Jl,2 I 


-n 


'32,31,2 


\m^—m • 


(96a) 

(96b) 


Again the nucleon-hole case differs from the above only for the contact 
terms, because a direct evaluation shows that the various £2 are just the 
limits of the £2* for M —> m. Only one contact term exists for the case , 
with Afg = —2/q^. 
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8 The spin 3/2 resonances 

We consider now the excitation of a nucleon to a spin 3/2 resonance (specif¬ 
ically the A(1232) ), assumed to be stable. The resonance is described by a 
vector-spinor held ip^ obeying the Rarita-Schwinger equations 


i (pS- 

= 0 

(97a) 


= 0 

(97b) 

dY, 

= 0 

(97c) 


(the last line is, more properly, a constraint) which can be deduced from the 
Lagrangian HH] 


£ = ^^|(i ^ + iuiYd'^ + d>^Y) 

+ ^ (3u; 2 + 2a; + 1) 7 ^ ^ 7 " + {Suj^ + 3a; + l) , (98) 


uj ( 7 ^ —1/2) being a free parameter. Each value of a; leads to the equations 
of motion (EZD, but does not prevent the occurrence of a spin 1/2 compo¬ 
nent in the a;-dependent vector-spinor 'ip^. Thus, to rule out these unwanted 
components, one usually introduces the projection operator on the spin 3/2 
space, which reads (in momentum space) 

R 2 = s"" - - A . (99) 

or, sometimes, its on-shell reduction 

_, pM^^ I 

-^3/2 -^3/2 lp2^M2 



The most common choices are to = —1/3, that leads to the Rarit a-Schwinger 
result, and a; = — 1, that corresponds to the Lagrangian 


£ = ^P, . ( 101 ) 

Another Lagrangian has been recently proposed, namely na 

£ = 7{-e'"'VTiap-Mg'‘7i^„ (102) 

in order to solve the so-called Velo-Zwanziger disease [iHiiin]. However we do 
not discuss such a Lagrangian here, as it describes a resonance propagating 


in an external electro-magnetic field with the (minimally coupled) AA 7 ver¬ 
tex, while we only consider non-minimal N-A transitions. Furthermore, the 
proposal (cni is seriously plagued by the occurrence of a pole at = M^/4, 
as the evaluation of the propagator (the inverse of F^*^) shows. 

Sticking to the more sound form (Ih 8 j) . observe that different values of 
uj not only alter the mixing between 3/2 and 1/2 spin, but also affect the 
off-shell behaviour of the A-hole propagator, that in fact reads 


p2 _ Y 
1 + UJ 

^ 6 M 2(1 + 2 a ;)2 


3' ' 




7V) 




3M"' " ' " ' 3M2^ 

|(1 + 3w)(7V + 7 V) + (1 -t- uj){Yp" - 7V) 


- 2MujYY + (1 + cu) /y V \ (103) 


Remarkably the choice a; = — 1 cancels in (uni all the terms with no analytic 
structure, which would otherwise contribute to the contact terms in the A- 
hole LF. 


8.1 The 0-index A-hole Lindhard functions 

Here the vertices have the form or where could be taken 

from Eq. (IMl) plus the non-conserved and Clearly 

fxyip, q) = Tr(/ + m)j^^Af^''jy^ (104) 

or, alternatively, 

fxyip, q) = Tr(/ + ■ (105) 

However, the vector 7 ^ is constrained by (inib) so that this current, when 
contracted with A^'^, is vanishing on shell, does not develop an imaginary 
part and consequently it can only contribute to the contact terms. The same 
occurs for 7 ^ 75 . Furthermore can be replaced by because 


Thus only the two currents 


JiiNAr = q^ 


(106a) 
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j2{NAr 


(106b) 


actually matter, at least for the part carrying analytic structure. Since we 
are dealing with a 0-index LF, the structure is given by (1771) . With self- 
explanatory notations the non-vanishing are found to be 


Concerning the contact terms, we will not give a detailed list for all the 24 
currents, because they all explicitly depend upon uj and display a double 
pole at a; = —1/2: hence they can diverge and the real part of the LFs 
becomes unpredictable. As a consequence, the RPA series based on the A- 
hole excitation (and, similarly, any calculation beyond the bare Free Fermi 
Gas) becomes unreliable. 

The same happens if we use the expression (Unsi) taking however the 
projection operator in the form (jlOOj) . We get indeed Eqs. (unzi) for the 
functions , but again the contact terms display a double pole in uj. 

Finally we can take the projection operator in the form TO . Now the 
expressions (unzi) are again valid, but the contact terms are independent of uj 
and, furthermore, they do not change in replacing with They dis¬ 

play however a factor {p + q)~‘^ coming from the projection operator Eq. (Ib9|l . 
For instance, in the current ji{NA)^ they take the form 



^ [ol 

Here the hrst term is just the function evaluated at M = 0, thus ex¬ 


plaining the introduction of the factor B^y{q'^)T^^\M = 0; |q|, go) in EQ- (EH). 
Explicitly the case (j/j/) requires 



(rn? — q^Y(m? -|- 2mM — q^) 


and furthermore Afo = {rri^ — 2mM — 2M^ -|- 3g^)/3 and Afo = 


(108) 
1/3. Finally, 









8.2 The 1-index function 


Now we consider the case of the transition from a scalar to a vector term, 
which leads to a 1-index LF. Here, besides the vectors discussed in the pre¬ 
vious subsection, we also need a tensor operator yielding a vector when con¬ 
tracted with the A propagator. Again one can set up a vast amount of 
tensors: we limit ourselves to those which give LFs with the same analytic 
structure, ignoring extra contact terms. 

We are thus left with eight possible currents, namely 

iv = + h-c. (110) 

with 


-pfiu 

^ M(V) 

2 Q+ 

(111a) 

■pflU 

^ E(V) 

= 


(iiib) 

^C(V) 

.3 fi + 1 
'2Q+Q_ 

2q^ (g V - P ■ (1(f) 75 

(lllc) 

^ M(A) 

.3 fi — 1 
'2 Q_ 


(llld) 

^ E{A) 

.3 /i- 1 
'2Q+Q_ 


(llle) 

^C(A) 

- ^2Q+q p-m'") 

(lllf) 

^S{V) 



(lllg) 

^S(A) 

= 9 ^^ 


(lllh) 


(here /i = M/m). We have followed in the above the work of Devenish et 
ah 1221 : these author show that, for a transition from a nucleon to a higher 
spin resonance (not necessarily 3/2), only three conserved currents enter the 
parity-conserving sector and as many the parity-violating one. In the low 
momentum regime tensors associated to these six currents correspond to the 
multipoles Ml, E2, C2, M2, El, Cl (the hrst three parity conserving, the 
other parity violating). We have added two other currents, which are not 
conserved, thus exhausting all the possibilities. 

We see that again LFs with only one vector index exists, but they only 
occur between a Coulomb multipole and the non-conserved currents (ITO . 
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The case {j 2 ^c(v)) contains only a in the integrand, thus the expression 
applies (hrst term in (IHTHl i with 


V 

42^rc(y) 

_ A O 

^ mM"^ 



. Mt (m2 

— q^y{w? — 2mM — 

32^C(V) 

" mM‘^ 

Q+Q- 

^10 

iMTq‘^ rri^ 

+ 2mM - 2M2 + 3g2 

m 

Q+Q- 

^20 

iMxq^ 


mQ+Q- 



,2\^2 


(112a) 

(112b) 

(112c) 

(112d) 


The case {j 2 ^s{v)) requires the full expression and we hnd (since = 

0 ) 


:i2^s(v) 

^10 






j 2 ^S(V) 

:i 2 '^S(V) 


'^20 


{q^ + M^- m^) 

{m^ — 2mM — — (f) 


{^2 , / 14-2 

3M2 
2 

3ii^ 

2 

3 

3M2 
1 

~3M2g2 

1 

^ ■ 


2^2 


— 2mM — q^){m^ — q^) 


The LFs with opposite parity simply obtain as 


ndi)M 

3iA x(A) 


(|qUo) 


ndi)M 

J2,i x(V) 


(|qUo) 



(113a) 

(113b) 

(113c) 

(113d) 

(113e) 

(113f) 

(114) 


8.3 The 2-indices functions 

8.3.1 Parity conserving-parity conserving Lindhard fnnctions 

We consider hrst of all the couplings Ml, E2 and C2 (parity conserving) 
in both vertices. These currents being conserved, we can directly apply 
Eqs. (gSl) and (HED- Hence the functions A^y and the corresponding 
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^M{V) 

^EiV) 

^C(V) 

3M|.g^ 

9M^g" 

3M^g^ 

m?Q+ 

m?Q+ 

7V?M^Q+ 



^E(V) 

^C(V) 

q 

9M# 

0 

Am?q^ 

4^2 g2 


Table 6: The A^{q‘^) functions for 
the vector-vector A-N parity con¬ 
serving currents 


Table 7: The functions for 

the vector-vector A-N parity con¬ 
serving currents 


A’s are needed. Remarkably, A^y and A^y are diagonal with respect to the 
channel indices and are quoted in Tables El and 0 respectively. Instead, a 
contact term coupling the multipoles Ml and E2 exists (it is reported in 
Appendix ini together with all the other contact terms). 

Both the multipoles Ml and E2 display the structure of Eqs. (gg) and iSD 
with, however, = 0, while in the Coulomb multipole C2 the contribution 
proportional to 11*'^ is absent. is given in Table 0 and furthermore 

Bay)C{V){(f) = - 2mM - q^) . (115) 

Observe also that in this sector flx^ 

8.3.2 Parity conserving-parity violating Lindhard functions 

This kind of LF is antisymmetric in the indices p, v and takes the form 

(see Eq. dSl). We thus need to specify the functions A^y and the contact 
terms. The only non-vanishing A^y are 

^E{V)M{A)i^'^) = ^^^{V)E(A)i.^^) = ^2^ 

and for the relevant contact terms (not A]^ = 0) we refer the reader to 

Appendix O Moreover, 


T-rO/i.h' _ T-rO/i.h' 

^^x{A)y{V) - ~'-^x{V)y{A) 


(117) 


(with x,y = M,E,C). 
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8.3.3 Parity violating-parity violating Lindhard functions 

The functions and follow from the parity conserving-parity 

conserving case with the replacements 


A 


E(A)E(A) 




- ^M(V)M(V)(^^) ’ 


•> 


A 


LT 

C{A)C{A) 


-A 


LT 

c(y)C{v) 





LT 

C{A)C{A) 


^C{V)C{V) \m^-m 


(118) 


As for the parity conserving-parity conserving case the rule (involving contact 
terms) holds. The contact terms (see Appendix 0 have a 

quite involved structure. 


8.3.4 Lindhard functions involving non-conserved currents 

A LF having the hrst vertex M{V) or E{V) and the second S{V) (this 
last corresponds to a non-conserved current) is non-vanishing and has the 
structure of Eq. iOD (a symmetric LF obeying the current conservation law). 
Thus it can be expressed in terms of Eqs. (P|) with 


o aS — aS — 

^^M{v)S{v) - ^E{V)S(V) - 2mq^ 


and 


q aLT — 

^^M{V)S(V) — ^E{V)S{V) 


QiMxq^ 


mQ+ 


The same happens for 11 


E{A)S{A) 


with 


A 


s _ 

E{A)S(A) - 


.5mQ+ 

*-^ 

mq^ 


4 LT _ 

MA)SiA) - , 


(119) 


( 120 ) 


( 121 ) 


while only displays contact terms (see Appendix E]). 

The LF has a different structure: it contains a symmetric 

current-conserving term, as in Eq. (gni), with = 0 uud 


A 


LT 

c{v)S(y) 


2iMTq^ 

mM‘^Q+ 






(122a) 
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B 


LT 

c(v)S(y) 


2iMTq^ 

mM^^Q+Q. 


q^){m^ — 2mM — q^) 


(122b) 


plus the following term proportional to q^ 




mM'^ 

^ iMj' 

^ mM^Q+Q_ 



q^Y{m^ — 2mM — q^)0*J {M = 0) 

+ '?o) 

m 


(123) 


with the contact terms given in Appendix O The symmetry relation 


_ _TT*(2)/21^ I 

^^C{A)SiA) ~ ■^^C(y)S(V) I m^—m 


(124) 


holds. 

The functions with initial vertex M{A) and E{A) and hnal vertex S{V), 
as well as '^*m^v)S(A) n 5 (y^s(yi) antisymmetric and are given by the 
second term of Eq. dH, namely 


with 


and 


qxm'^g^pe^^^^U:l_{M;\q\,qo) , 


^M{A)S(V) 

^M(V)S{A) 


'^^E{A)S{V) 

2Mj' 



m 


A 


S(V)S(A) 



(126) 

(126) 


(127) 


while ^eIv^)s{A) contact terms and = 0. 

Next we consider the symmetric function Since the current 

is not conserved, it is contributed to by all the terms in Eq. (jlTj). Thus it 
displays a conserved part, which has 


A 


s 

SiV)SiV) 


4Q_ 

3g2 


(128a) 
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aLT 

^S(V)S{V) 


(128b) 


tjLT 

^S{V)S{V) 


4Q_ 

3M2 

4(q'^ — + 2mM) 

3iW^ 


(128c) 


plus a non-conserved one (last two terms in (jUj)). The latter requires the 
knowledge of - in turn fixed by (see Eq. (jSH)) 


^S{V)S{V) 

TjV 

^s(y)s(v) 


2Q_ 

3M2g2 

2 (m^ - 


(M^ — + q^) 

- q^){w? — 2mM — q^) 


3M^q^ 


(129) 

(130) 


Finally for the third term, which, being proportional to q^q’^, is associated 
to a scalar quantity, as in Eq. (EH), we get 


~aS _ V-V+ 

^5(y)S(n) - 

im? — q^Y{m? — 2mM — q^) 

Bs(v)S(v) = 

Again the symmetry relation 

tt*(2)aii^ _ _ Yr*{2)^u I 

^^s{A)s{A) ~ ■^^5(y)S(y) I m—m 


(131) 

(132) 


(133) 


holds. 

Note that the interchange of the initial and final vertex entails an inter¬ 
change also of the indices p and z/. 


A The elementary functions 

The functions defined in Eq. (ESI are given by 


ry-i*[l,0] 

^5+ 

= 2Tt°' 

(134) 

r^*[2,0] 

^5+ 

= 2gV2^'°’ 

(135) 

r^*[2,0] 

^S- 

= 

(136) 

rTi*[3,0] 

= -H(4m^|q|^ - 3y,‘)Tl"l + 5(3,„^ + |q|^)SPl 

(137) 
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(138) 

(139) 


rT-i*[4,0] 

^5+ 

Tn*[5,0] 

^5+ 


ry-i*[l,0] 

^V+ 

rT-i*[ 2 , 0 ] 

^V+ 

rT-i*[ 2 , 0 ] 

^V- 

rT-i*[3,0] 

^V- 

rT-i*[4,0] 

^V- 


ry-i*[l,0] 

^L+ 

ry-i*[ 2 , 0 ] 

^L+ 

rT-i*[3,0] 

^L+ 


^T+ 

ry-i*[ 2 , 0 ] 

^T+ 


2pq^{p^q^ — 4m^|q|^)Tf°' + 8pq‘^{3ql + |q|^)Tf^] 


/ 4 2 3 2| |2\2 64 4 4 

[q p —Am q )- m q 

5 




64 


16gV^(3g^ + IqD - y^^|q|^(5go + IqD 




32 


+ IqH 


2 4, 


|q| 




|q| 




-2p|q|''rw 

(m^T^ - 4X12]) 

2 p|q|'T*'!’°' 

Am‘^qo\qf 


q'^p^ -m^|q|^ ) 

5 


16go|q| 


2 r 


( 9o+glql 


2 / „2 


7, 




32go|qr 


^o + s'qi 


T[4] 


+4|q|")XP'] 

6q 


2 r7n*[l,0] 

<1 pTl+ 


2 ^ 2|^|2 r 


2 m^q| 
■3 3,4 


Q 


2 |q| 


3 q^ 
8 |q| 


A 9 72 91 12 

q p -m |q| 

5 


2 :[o] 


i2 r 


(9? + 3 |qr)<,V - + 5 |qr)(|qr + 59S 

9‘ + f 9„"9"] X|4] 


3 q^ 

_1^22;[0] ^ lij[2] 

3 3 

fpTpP' 


(140) 

(141) 

(142) 

(143) 

(144) 


(145) 

(146) 

(147) 

2 :[ 2 ] 
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(148) 

(149) 



(150) 


ry-i*[3,0] 

-^T+ 



4 2 ^ 21 

q p — -m q 
5 


2 


Xio] 


+ 

+ 


4 

3 

16 

y 


gV -4m^ ( + -|q 


X [21 


-lo + jlql 




(151) 


B The elementary functions 


r_ jjtfoi A *rii 

In Sec-IHwe have derived the functions . Actually we also need and 
*[2l 

, which be be provided in this appendix. It is clear from the derivation 
of that the entire class of functions defined by Eq. (ITH|l has a common 
structure, namely 

= a(r' + + a'^'k + affe • (152) 


The coefficients are given by 


[0] 

«0 

= 0 


[0] 

a[ 

P 


IGtt^ 

[0] 

«2 

_ qop + 2Ef 

327r2 q 

[0] 

«3 

A* 


647r^ 

a™ 

pkp 

327r2 


qo 

327r2 


«[■! 

1 


1287r2|q| 

4" 

qo 

pA* 

1287r2 

4=’ 

_ kpqo 
487r^ 

P^- 


2rrr 

g 2 

44-(|q|' + %V + 


Am^qQ 


pEp 2m? 
2 go 




(153a) 

(153b) 

(153c) 

(153d) 

(153e) 

(153f) 

(153g) 

(153h) 

(153i) 
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a 


a 


[ 2 ] 

1 

[ 2 ] 

2 



P 


1927r2 

1 

3847r2|q| 


/(|q|' + 3go^) 


6 m^ 


ihf + Ql) 


PQo 


(3|qr + goV- 


2. 2 


+ 8E'j, 


A* 4m^|q|^ 


7687r2 


-(|qr + 3goV 


(153j) 

(153k) 

(1531) 


C The contact terms for the 2-indices A-N 
functions 

In this Appendix we list the 69 non-vanishing contact terms associated to 
the two-indices LFs of Sec. 13.41 We use the notation ^mnXi{Yi)X 2 (Y 2 p 'where 
a = S', V, LT is associated to the Lorentz components, Xj = M, E, C, S to 
the magnetic, electric. Coulomb or scalar nature of the two currents {i = 1, 2 ) 
and Yi = V,A to the their vector or axial parts, respectively. They are given 
by: 


1) vertices M(y)M(y) 

3M| / 

AoMivmv) - ^ J 

■^20M(y)M(v) = “4^2^2g^ (Q- ~ ~ 

xS 3M|. 

AoM(V)MiV) - -4^2^2g^ 

\ LT 3M^q‘^ 

AoM{V)M(V) = ^2Q2 

2) vertices M{V)E{V) 

3M|. l' AnY'^ 

^mM(V,EiV> = h “ 

3M^ 

>'loM(v}E(r) = 4n,2^2Q2Q_ ( 39 V^ - 4 mM()V - imHlrSm) 


(154) 

(155) 

(156) 

(157) 

(158) 

(159) 
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■^30M{V)E{V) 


[3(5+ — 4m(2m + M)] 


(160) 


A 


s 

40M{V)E{V) 


-^10M{V)E{V) 


X 


LT 

20M{V)E{V) 


3M| 

Am?q^Q\Q- 

3M| 

Aim?q^Q\Q_ 

JMMLin 

rn^QlQ^ 

3M^q‘^ 

m^QlQ. 


Arn^) 


(161) 

(162) 

(163) 


3) vertices E(y)E(y) 


■^10E{V)E{V) 

■^20E{V)E{V) 

■^30E(V)E{V) 

■^40E{V)E{V) 

■^50E(V)E{V) 

■^10E{V)E{V) 

■^20E{V)E{V) 

■^30E{V)E{V) 


3M^ 

(164) 

imWQlQ- ^ 

20rn^ MTSm\ 


(165) 


(166) 

4m(13m — 3M)q^p — 4m^(9m + 7M)5m] 


3M^ 

(167) 

3M| 

im?q^Q\Q‘^_ 

(168) 


(169) 


(170) 

12M^q‘^ 

m^QlQl 

(171) 


4) vertices C(V)C(V) 


X 


LT 

ioc(y)c'(y) 


A 


LT 

2oc(y)c(y) 


lMhL(n 

m^QlQl 

3M^q^ 

~w?QlQl 


3M^ + 4g2) 


(172) 

(173) 
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5) vertices M(y)M{A) 


y 3iMT6m 2 

^20M{V)M(A) = “^2g2g_? P 

^ y 3iMTSm 

hoMiV)MiA) = 2m^QlQ_ 


6) vertices M(y)E{A) 

^ y 3iMT5m 

^20M(V)EiA) = -2m^QlQ_ 

^ y 3iMjAm 

^30M(V)E(A) = -2^2q2^q_ 


[3(5- + 4m(2M — m) 


7) vertices E(y)M{A) 

^ y QiMTSm 

^20E(y)M(A) = ^2Q2 ^q 2_ 

X {—+ m{5m + 3M)q^p + 2m^g^} 

^v 3iMT5m _ , _ 

hoE{v)M{A) = ~2 w?Q\Q'^ 4m(2m + M)] 

3iMTSm 
\oEiv)MiA) = 2m^QlQl 


8) vertices E(y)E{A) 

>^20Eiv)E{A) = -2 ^^/g2 [3gV - 4m(5mgV - Am\^] 

^30E{V)E{A) = ~ 2m?QXQ"^ {3Q+ - 4m(2m + M)} 

^ y 3iMTSm 

\oEiv)EiA) = -2m^QlQl 


9) vertices M{A)M{A) 

3(5TTT-^ 

•'foM(/l)M(/l) = ~ 4 „ 2 ^ 2 |g_ P' 3 - - - 3 V)| 


(174) 

(175) 

(176) 

(177) 

(178) 

(179) 

(180) 

(181) 

(182) 

(183) 

(184) 
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>^WM{A)M{A) = 2mWQ+Q^ + 4m(3m + M)q^p + Arn^q^] 

?>5w? 

30M(A)MiA) = ~2m^q^QlQl 

X {3g^p^ — 4m(2m + M)q^p + Am?MT^^m + M)} 

^%M{A)M{A) = ~4m‘^q2Q2^Q2_ {^<5+ “ 8mMr} 

3(5m^ 


>^loMiA)MiA) Am^QlQl 

^LT ?>5m?q^ 

K0MiA)MiA) = -^2g2g2 

X {3g^p^ — 4m(4m + 3M)q^p + Am^MTihm + 3M)} 
35m?q^ 


'^20M{A)M{A) 


A 


LT 


m2g2g2_ 

35m^q^ 


{g_ - 4m^} 


30M(A)M(A) ^2g2g2_ 


(185) 

(186) 

(187) 

(188) 

(189) 

(190) 

(191) 

(192) 

(193) 


10) vertices M{A)E{A) 


mOM{A)E{A) 


'^20M(A)E{A) 


'^30M{A)E{A) 


\S 

MoM{A)E{A) 


■^50M{A)E{A) 


-^10M{A)E{A) 


'^20M(A)E{A) 


3(5m^ 

4m2g2 

35m? p 
m^Q+Q- 


35m? 

2m^QlQ? 

35m? p 
~w?QlQ? 
35m‘^ 

Am^q^QlQ 

35m?q^ 

m^Q+Q- 

&5m?q'^p 

~m^QlQl 


{3gV 


2 


2 


4m^} 


(194) 

(195) 

(196) 

(197) 

(198) 

(199) 

( 200 ) 
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A 


LT 

30M(A)E{A) 


36m‘^q‘^ 


( 201 ) 


11) vertices E{A)E{A) 


-^10E{A)E{A) 

■^20E{A)E{A) 

-^30E{A)E{A) 

-^AOE{A)E{A) 

■^50E{A)E{A) 

-^10E{A)E{A) 

-^20E{A)E{A) 

■^30E{A)E{A) 


35w? 


Am?q^Q- 

35m? 


2m^q^Qj^Q‘? 


[3Q- + 3m5m] 

[Sg^ — 4m (2m — M)q^p 


36m^ 


2m‘^q^Q\Q'? 

35m? 

Am‘^q^Q\Q? 

35m? 

Am^q^Q?Q\ 
35m?q^ 


[5gV^ — 4m (2m + M)q^p 


[bQ- — 4m(2m — 3M)] 


m?Q+Q? 

35m?q^ 

m^Q?Q\ 

35m?q^ 

'm?Q?Ql 


[SQ- + 4m(5m] 
[3(5+ — 4mMT] 


( 202 ) 

Am? MT5m\ 

(203) 

4m^M7’(5m] 

(204) 

(205) 

(206) 

(207) 

(208) 
(209) 


12) vertices C{A)C{A) 

^WC{A)C{A) = ^ 2 Q 2 q2^ [ 3 ( 5 + - (2m + M)^] 

35m? q'^ 

A20C(A)C(A)^2g2_g^ 

13) vertices M{V)S{V) 

xs iMrp 

5iOMiV)SiV)- ^ 

^ g iMj^ 

^20Miy)SiV) = 


( 210 ) 

( 211 ) 


( 212 ) 

(213) 
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14) vertices E(y)S(y) 


^io£;(y)s(v) 


A 


iMj 


mq^Q+ 

iMj' 


A 


20E(v)siv) - 2mq^Q+Q. 

= mq^Q+Q. 

^lt dtiMxq^ 

KoEiV)SiV) = 

15) vertices C{y)S{y) 

2iMTq^ 


[Q+ - QraMr] 

[IIQ- + I6m5m] 


^ioc'(y)s(y) 


^2oc'(y)s(y) 


mQ+Q_ 

iMj' 


mQ+Q_ 

16) vertices M{A)S{V) 

^20MiA)SiV) - 

17) vertices E{A)S{V) 


-^20E{A)S{V) 


6m 


mQ+Q- 
18) vertices S{V)S{V) 

4 


X-^ 

^ios(y)s(y) 


x*^ 

^ios(y)s(y) 


X*^ 

^ios(y)s(y) 


3g2 

+ 2mM — 2M^ + 3g^ 


3g4 


1 


(214) 

(215) 

(216) 
(217) 


(218) 

(219) 


( 220 ) 


( 221 ) 


( 222 ) 

(223) 

(224) 


44 



19) vertices E(y)S{A) 


_ 2Mt 

^20E(V)S(A) - ^Q^Q_ 

20) vertices M{A)S{A) 

\S ■ 

-^10M(A)5(A) = 

2i5mp 

^20MiA)SiA) = -^Q^Q_ 

^30M(A)S{A) = ^q2Q^Q_ 

Ai5mq^ 

Aiom(A)5(A) = 

21) vertices E{A)S{A) 

StDj 

AfoB(A)S(A) = 

AfoB(A)s(A) = -^ ^q 2 Q^Q_ [3<5- + 4m(2M - m)] 

\S 

^30EiA)SiA) = -^p^q2Q^Q_ 

.LT _ 

^10E(A)S(A)- 


(225) 

(226) 

(227) 

(228) 

(229) 

(230) 

(231) 

(232) 

(233) 
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0 0.5 1 1.5 

q(GeV/c) 


Figure 1: The singularities and of the function U*li in the (Iqlj^o) 
plane. Here, having in mind nucleons and A’s, we have m = 0.938 GeV and 
M = 1.232 GeV. Solid lines denote the four singularities Q±± and dashed 
lines the ±Q±. Dotted lines denote the light cone, kp = 0.264 GeV/c. Light 
grey (yellow if in colors) region: Q^_ > qo > medium grey (cyan) 

region: —Q^ > go > Q^_', dark grey (green) region: > Qq > very 

dark (red) region (not visible in this hgure): > qo> Q++. 
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Figure 2: As fig. ^ but with kp = 1.5 GeV/c, in order to emphasize the 
region > go > Q++) otherwise invisible. 
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